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Abstract 

I derive a single factorization formula which summarizes all soft and collinear diver- 
gences of a tree-level gauge theory amplitude. The singular factor squared is in a certain 
sense the generalization of the Catani-Seymour dipole factorization formula. 
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1. Introduction 



The factorization properties of gauge-theory amplitudes in infrared-singular limits have played 
an important role not only in our conceptual understanding of infrared safety, but also in explicit 
calculations at next-to-leading order in perturbative QCD. These infrared singularities may be 
classified into two types: soft and collinear. The former arises when a gluon four-momentum 
vanishes, k s — ► 0; the latter when the momenta of two massless particles become proportional, 
k a z(k a + kb), kb — ► (1 — z){k a + kb). The properties of non-Abelian gauge-theory amplitudes 
in these limits are easiest to understand in the context of a color decomposition [1]. For tree-level 
all- gluon amplitudes in an SU(N) gauge theory the color decomposition has the form, 

A t ™ e ({k i ,\ i ,a i }) = £ Tr(T^) • • • T a ^)A^(a(l^ , . . . , n A ")) , (1.1) 

where S n /Z n is the group of non-cyclic permutations on n symbols, and j Xj denotes the j-th 
momentum and helicity. As is by now standard, I use the normalization Ti(T a T b ) = 5 ab . One can 
write analogous formulae for amplitudes with quark- antiquark pairs or uncolored external lines. The 
color-ordered or partial amplitude A n is gauge invariant, and has simple factorization properties in 
both the soft and collinear limits, 

AT e (- . . , a\ &\ . . .) Yl SpUt t _T(a A «, b x ») A* e _\(. . . , (a + b)\ . . .) , 

A=± (1.2) 

• • , a, s\ b, . . .) Az^ Soft tree (a, s x ° , b) A^ e _\ (. . . , a, b, . . .) . 

The collinear splitting amplitude Split tree , squared and summed over helicities, gives the usual 
unpolarized Altarelli-Parisi splitting function [2]. While the complete amplitude also factorizes in 
the collinear limit, the same is not true of the soft limit; the eikonal factors Soft ree get tangled 
up with the color structure. It is for this reason that the color decomposition is useful. Its use, 
and the simple factorization properties of the partial amplitudes, allowed Giele and Glover, in their 
pioneering paper [3], to derive simple and universal functions summarizing the integration over soft 
and collinear radiation. Of course, it is possible, as was done by later authors [4,5], to square the 
amplitude first, and tangle up the color factors with the eikonal factors, thereby recovering earlier 
forms of these two factorizations [2,6], but this obscures the clean structure in eqn. (1.2). 

The color decomposition (1.1) is also the form in which one writes tree-level amplitudes in an 
open string theory, and the gauge-theory partial amplitudes may be computed as the infinite-tension 
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limit of the corresponding string amplitudes. These are given by the Koba-Nielsen formula [7], 

„ n—3 

^KN = T n-4 2 -„/2 j Yl d X% (x b - X a )(x n _ 2 - X a )(x n _ 2 - X b ) 



7 = 1 

X a <X 1 <X b ---<X n _ 2 

(L3) 

x I | \xj - Xi\ T ki ' kj/2 exp 

i^jZzE j multilinear 



El £j • £j T~ki ■ £j 



2 (Xj X^) (Xj Xj 



which I have written in a form that will be useful for the derivation to be presented in this paper. 
In this equation, {ki,£i} are the momenta and polarization vectors of the gluons; E denotes the 
set of all external legs, {a, 1, 6, 2, . . . , n — 2}; r is the square root of the inverse string tension; and 
the subscript 'multilinear' tells us to extract the terms in the exponential linear in each of the 
polarization vectors. Three of the Xj may be fixed at will; it will be convenient to choose x a = 
and Xb = 1. The gauge theory limit is given by r — > 0. This formula gives an explicit form for the 
n-gluon amplitude, and examining its limits is thus a convenient method of deriving [8] factorization 
formulas such as those in eqn. (1.2). 

More recently, Catani and Seymour [5] wrote down a dipole factorization formula which gives 
a single function capturing the singular behavior of the squared matrix element in both the soft 
and collinear limits. Of course, since it is at the level of the amplitude squared rather than the 
amplitude, and since it does not take advantage of the color decomposition, it is not quite a true 
factorization but is still tangled up with the color algebra. The purpose of this note is to derive 
a single function unifying the soft and collinear factors in eqn. (1.2) at the level of color-ordered 
amplitudes. This will provide a true factorization. I will discuss its connection with the Catani- 
Seymour dipole factorization formula later. The derivations in ref. [8] treat the soft and collinear 
limits independently. 



2. Derivation of a Unifying Factorization Formula 

Before embarking upon the derivation, I note that the Koba-Nielsen form is valid whatever 
the dimensions of the external momenta and polarization vectors, so long as the former are on- 
shell, and the latter transverse. It thus allows for a treatment corresponding to any of the possible 
variants of dimensional regularization: the conventional scheme (CDR) [9], the original 't Hooft- 
Veltman scheme (HV) [10], dimensional reduction (DR) [11], or the four-dimensional helicity scheme 
(FDH) [12]. 

The forms of the collinear and soft factorization functions suggest that we can associate the 
singular limits with a trio (a, 1,6) of color-ordered partons: either a soft parton 1 with two hard 
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neighbors (a, 6), or a pair of collinear partons neighboring another hard parton, a || 1 next to b or 
1 || b next to a. This in turn suggests generalizing the derivation in ref. [8] to extract all singular 
behavior in either of the two invariants s a i or su,. This will summarize the singular behavior 
whenever the trio of momenta become degenerate without either a or b becoming soft, that is as 
A(a, l,b)/s^ b — > 0. (A(a, 1,6) is the Gram determinant of the (a, 1,6) system.) In contrast, the 
derivation of the soft limit in ref. [8] takes both s a i and su, to be small, and the derivation of 
the collinear limit takes only one invariant (say s a i) to be small. The derivation below otherwise 
follows the same general lines: study the singular behavior of an amplitude, and identify the parts 
of the singular expansion that give rise to an amplitude with fewer external legs. 

Let us then evaluate equation (1.3), focussing first on the singular behavior in either of the two 
invariants s a \ or su,. Within this singular limit, we may also simplify expressions by expanding 
around the field-theory limit r — > 0. If we examine the equation, we see that singularities in these 
two invariants can arise only from the integration region x\ ~ j;„ = or ij ~ Xf, = 1, and then 
only in those terms which have a lone inverse power of X\ or 1 — x\ coming from the expansion 
of the exponential. (The integrals of x~ 2+s or (1 — x)~ 2+s are finite by analytic continuation as 
5 — > 0.) Separating out all terms containing the integration variable X\, we obtain 



r n-4 2~n/2 



I 



n — 3 



dXi (x b - x a )(x n - 2 - x a )(x n - 2 ~ X b ) 



X a =0<Xl<X b = l---<X„-2 



1=1 



X X 



X < 



t k a -k\ 



(i - Xl y k "- ki Yl fa - xtf ki - ki 



cxp 



+ 



jG{2,...,n-2} 



rk a • £i rfcj, • £i rfei • e a rfei • e b 



(1-Xi) 2 xi (1 — xi (l-^i) 



+ E 



ie{2,...,n-2} 



£j • £i 



+ 



(Xi-Xl) 2 (Xi-Xl) (Xj-Xl) 



r. - Xi\ T2ki ' kj/2 exp 



El £i ■ £ 
9 



2 ^i) (^i Xj 



multilinear 



(2.1) 

where H = {a, 6, 2, . . . , n — 2} denotes the collection of external legs excluding leg 1. The terms 
on the last line form the integrand of the (n — l)-point amplitude, with leg 1 excluded. These 
will carry along with them (n — 5) powers of r, leaving one to be associated with the x\ integral. 
Noting that we can partial-fraction terms containing inverse powers of both x\ and (1 — xi), we 
see that the x\ integral itself can produce at most two inverse powers of r. (Each inverse power of 
an invariant comes along with two inverse powers of r, but the invariants cannot appear together 



in the denominator of any given term.) Thus terms in the integrand containing more than a single 
additional power of r will disappear in the gauge-theory limit. This disposes of the IT/e{2 n-2} 
factors on the second line, after expansion in powers of x\. Using this observation in expanding 
the first exponential as well, in addition to dropping terms that have no inverse powers of x\ or 
(1 — x\) or are not linear in £1, we find, 



T n-4 2~n/2 



ra— 3 



X a =0<X 1 <X b = l---<X„- 
( 



i=l 



X < 



i- E 



Y\ dXi (X b - X a )(x n - 2 ~ X a )( Xn -2 ~ X b ) X[ ka ' kl (l - Xl ) T ^ 
Tki ■ Si 



(xj-xi)J [ x\ (l-Xif 



ie{2,...,n-2} 

rk\ -£ b e a -£i rfci • e a £bj_£i _ rk a ■ gi rk b ■ £\ 



(1 — xi) x\ X\ (1 — X\) 2 X\ (1 — Xl) 



+ 



- n 



rki ■ e a rk\ ■ e b 



Xi (1-Xi)_ 



E 



£i ■ £i 



ie{2,...,n-2} 



(Xi - Xi) 2 



exp 



s.p. 



E 



ie{2,...,n-2} 



rki ■ £j 

(Xi - Xi) 



it" k 2 ' k j I 2 



exp 



E 



multilinear 



(2.2) 



The terms with two inverse powers of x\ or 1 — x\ will not give rise to poles after the standard 
analytic continuation in the invariants. We can thus drop them inside the first set of brackets. This 
subtraction is indicated by the subscript 's.p.'. (It is equivalent to freezing xi at either or 1 as 
appropriate, in the exponential on the penultimate line.) 



Partial- fraction to isolate the singularities in the X\ integral and expand the non-singular 
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factors as appropriate either around x\ ~ or x\ ~ 1, 



T n-4 2~n/2 



/ 



n — 3 



x a =0<a;i <£(, = !• ••<a;„_2 



i=l 









1 


< 









]]_ dxi (x b - x a )(x n - 2 - x a )(x n - 2 - x b ) x\ ka - kl (l - xi) T k »- kl 
rki ■ Si 



+ 



l-rkx-Sb- ^2 

ie{2,...,n-2} ~* 

1 + rk x ■ e a - ^2 



+ 



1 



rk-y ■ £i 



Xl 



rki -£b- y] 



iG{2,...,n-2} 

rki - Si 



E b ■ £\ 



ie{2,...,n-2} 



fo-l)J 
e a • £i + rki -£ a eb-£i - rk a ■ e\ 



Ei ■ £i 



X? 



+ 



i€{2,...,n-2} 

+Tki -£ a + ^ 

ie{2,...,n-2} 



rki ■ £j 
{xi - 1) ! 



-rki -Eh 



ie{2,...,n-2} 



Ei ■ £\ 
(Xi - If 



E b -£i - Tki ■E b £ a -£i+ Tk b ■ £\ 
Tki ■ Ei 



exp 



E 



ie{2,...,n-2} 



(Xi ~ Xl) 



- n 



xj - Xi\ T ki ' k i/ 2 exp 



^""^ 1 Ei ■ £j ^ Tki ' £j 



multilinear 



(2.3) 



It is also convenient at this stage to extract all the terms linear in £ a or £ b from the exponential 
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on the last line, whence we obtain 



n — 3 



r n ~ 4 2""/ 2 J 11 d Xi (x b - x a ){x n - 2 - x a )(x n _ 2 - x b ) r x\ k ^ (1 - Xl ) T kb ' kl 

X a =0<X\<Xb = l---<X n -2 

ijije{2,...,n-2] 



X < 



< 






Xi V 



x\{x 3 - If 



■ ei + k\ ■ 



)2 



+ 



(1-Xl) 



ie{2,...,n-2} 



*?(*, - 1)2 

k\ ' S-i E a * £j 
(Xi - l) 2 X 2 



+ h ■ e a ^2 



x*( Xj -l)\ 
e b • ei - k\ ■ e b e a ■ e\ + k b ■ E\ e a ■ e b 



+k b ■ ei ^2 



x exp 



- £ 

ie{2,...,n-2} 



X 2 ( Xj -l) 2 



rk\ ■ Ei rk a ■ Ej rk b ■ e, 



/,., ^ 



i£{2,...,n-2} 



- n 



\ Xj - Xi \ T ki ' kj/2 exp 



(Xj-Xi) Xi (Xj - 1) 



1 Si ' £j ^ • Sj 

_ 2} 2( Xj -Xi) 2 [Xi Xj 



multilinear 



(2.4) 



Next, introduce polarization vectors £a,£ b for two new momenta and replace two dot 

products in each term via the identity 



■«"" = - £ 

polarizations h 



£ (hf(h) 



(2.5) 



where the sum runs over both physical and unphysical polarizations. The values of k b will be 
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fixed later; for the moment, they are arbitrary. Collecting terms, we obtain 



r n-4 2~n/2 



/ 



71 — 3 



Y[ dXi (X b - X a )(x n - 2 ~ X a )(x n - 2 - X b ) T X[ ^ (1 - X X ) T fc *" fcl 



x a =0<xi<x b = l---<x„ 



i=l 



polarizations(d,o) 



^-£g ■ £6 ^ - £ a ■ £l fei ■ £a + &i • £ a £ a • £i ~ &a " ^1 " £a) 



g a ' £ i £ b' £ 3 



exp 



- E 

ie{2,...,n-2} 



xnx,- — 1) 

i/jG{2,...,n-2} * V J 7 
+ r- ] __ x \ £ a ■ £ a (^b " £l " £g ~ " £ fe £g • E\ + A; fe • £i £ b • £g) 



(Xi ~ 1) 



£ a ' £i £ b ' £ 3 
xf(Xj — l) 2 



exp 



- E 

ie{2,...,n-2} 



Tfcq • £j T(/C b + fcl) • £j 



(Xi - 1) 



>< n 



xj - Xi\ T ki - k il 2 exp 



E 



multilinear 



(2.6) 

where 'multilinear' is now understood to refer to £ and e* independently, and where I have set x\ 
in the exponential to or 1 as appropriate to the pole structure of each term. 



Take x- a = 0, x b = 1, and define H = {— d, —b, 2, ... ,n — 2}; '—a' signifies the replacement of 
(ka,E* & ) by (fc_a,£_a) = (—ka,e*~). We can absorb the factors in brackets into the exponential on 
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the last line, 



r n-4 2~n/2 



n — 3 



I J dx, (x 6 - X a )(x n _ 2 - X a )(x n _ 2 - X b ) T x\ k ^ (1 - X X ) T ^ 



x a =0<x 1 <x b = l---<x n _ 2 



1=1 



polarizations(a,6) 



■ £b [~ e a -£iki-£a + fa-£ a £ Si -£ 1 - k a ■ £\ £ a ■ £„) 

T~{k a + fa) • £j | Tk b -£j | Tk^-Ei | rfeg • £j 



x exp 



E 

iG{2,...,n-2} 



+ 



+ 



+ 



(Xj - 1) X; (Xj - 1) 

+ T^-^TT^a • £ a ( £ b ■ ?! h ■ £~ b ~ h ■ £ b £g ■ E\ + fc 6 • £i £fc ■ £g) 



x exp 



E 

iG{2,...,n-2} 



rk a ■ £j r(k b + fa)-£j rfcg • £j Tfcg • e 3 - 



- 1) 



(Xi - 1) 



>< n 



IT" /c^ ■ Alj / 2 



exp 



^""^ 1 £j • £j ^ T/Cj • £j 

2 (xj x^) (xf 

' •■ j ■ a 



multilinear 



(2.7) 



up to terms which vanish in the gauge-theory limit. 

If we now choose kg, = —k a — ki and k b = — k b in the first term, and ka = —k a and k b = —k b — ki 
in the second, we can perform the Xi . . . x n s integrations. They yield an (n — l)-point amplitude, 
so that the singular limit of the original amplitude becomes 
-2 r 1 



T 

V2 



f dx 1 x? ka - ki (i-x 1 y ak "- ki 

Jo 



*E 

polarizations (a, b) 



-^-£ b -£b ( - £a ■ £l fa ■ £ a + fa-£ a £ a -£ 1 - k a ■ £i £ a ■ £„) 



xA™ 1 {...,-a = a + k 1 ,-b = b,...) 



(2.8) 



+ JY^pj £ a ■ e a (eb ■ £l fa ■ £g - fa ■ £ b £g • £i + k b ■ £\ £ b ■ £g) 



x A^ 1 (...,-a = a,-b = b + fa,...) 



Performing the xi integration and taking the r — > limit then gives us 
1 



E 

polarizations(a,b) 



Sal 



£ b ■ £b V 3 (a, 1, a) A™\ (. . . , -a = a + fa, -b = b, . . .) 



-— £a-£ a F 3 (6,l,6X r f e 1 (...,-a = a,-6 = 6 + fc 1 ,...) 



(2.9) 
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where Vs(a,b,c) is a certain form of the color-ordered three-point vertex (up to a factor of i), 



V 3 (l, 2, 3) = \/2 (-£1 -e 2 k 2 -e 3 + k 2 - e x e 3 -e 2 - &i ■ e 2 e x ■ e 3 ) . (2.10) 

In the singular limit, one might write down other forms of the three-point vertex. However, later 
formulae will not conserve momentum amongst the three momentum in these vertices (they will 
of course conserve momentum overall), and as we will be taking the ratio of two quantities which 
both vanish in the singular limit, we must be careful which form we use. (With this in mind, one 
could also use the recurrence relations formalism [13] to derive eqn. (2.9).) 

If we can find a pair of reconstruction functions k~ b = f &b (k a ,ki,kb) such that {k a ,k^) — > 
— (k a + ki, kb) sufficiently quickly when k\ \\ k a , and likewise (k a , k b ) — > — (k a , kb + k\) sufficiently 
quickly when k\ \\ kb, then we can pull out a common factor of the surviving hard amplitude. 
Furthermore, gauge invariance allows us to replace the sum over polarizations with a sum over 
helicities, so we obtain for the singular limit, 



helicities(a,fo) 



—e b -e b V 3 {a,l,a) - — e a ■ e a V 3 (b, 1, b) 

Sal $bl 



A*l\(...,-a(a,l,b),-b(b,l,a),...). (2 n) 



(In the four-dimensional helicity scheme, only the ± helicities are needed; in the conventional and 
't Hooft-Veltman schemes, 'e' helicities are needed as well.) The factor in brackets is the antenna 
factorization function, 

Ant(a,6< — a, l,b) = — e, ■ e h V 3 {a, 1, a) - — e a ■ e a V 3 (b, 1, b) . (2.12) 

Sal Sbl 

This function factorizes the singular region of three-particle production onto the antenna 
formed by two hard particles. It depends on the helicities of the original particles, as well as 
the helicities of the resulting hard particles. 



3. Explicit Forms 

In the k a \\ k\ collinear limit, the second term in eqn. (2.12) is non-singular, and we recover 
the usual collinear splitting amplitude, given by the first term [8] (note that the normalization 
conventions there are slightly different from those in the present paper). Similarly, in the kb \\ &i 
limit, only the second term contributes. In the soft limit (ki — ► 0), /e a = — k a , k b = —kb, and we 
can neglect all terms proportional to k\ in V 3 ; we then obtain 



' ^-a ^6 ' ^b 



-—k a -e! + —k b -e 1 

Sal Sbl 



(3.1) 
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which is just the usual eikonal factor. Thus we recover known factorization formulae in the appro- 
priate subregions. 

Sensible reconstruction functions satisfy momentum conservation everywhere, and also keep 
k~ g massless everywhere (not just in the singular limit). This corresponds to taking the sum of 
the original momenta, K = k a + k\ + kb, and rewriting it as the sum of two massless momenta. 
Working in the frame where K = 0, we see that there is a degeneracy corresponding to the relative 
angle of b and — b. In the singular limit, this degeneracy is fixed by the requirement that the angle 
vanish; but away from the limit, it survives. This means that the choice of reconstruction function 
is not unique. 

One possibility is to use the Catani-Seymour choice, 

K 2 



k a — k a k\ + 



Sal 
Sib + Sab 



-h, 



h = - 



Sib + S a b 



(3.2) 



which however is asymmetric in our case (and is not valid for kb \\ ki). A more symmetric choice 
is to take 



s bb = ~ 



SglSlb 
Sib + S ab 

SgjSlb 
Sla + Sab 



n 



2(1 -n) 
P+l 



(1-ri; 



1 - 



2ri 
P+lJ ' 



(3.3) 



where 



Sib 
Sla + Sib 



4n(l - ri)si a s u 



'1 + 



K 2 S a b 



(3.4) 



(Note that p — > 1 in all singular limits.) This is equivalent to 



k a — 



1 + P + 
l-p- 



sibjl + p- 2ri) 

Sla + S ab 

£ib(l + P~ 2ri) 

Sla + S ab 



1-P + 



Sla(l 



p-lTi) 



k a 



{i-n)ki-- 



l + P 



h 

Sib + Sab 

sia(l ~ P~ 2ri) 

Sib + Sab 



(3.5) 



kb 



These approach the appropriate Catani-Seymour choice in both collinear limits, a \\ 1 and b \\ 1. 

So long as we stay in four dimensions, we can use the spinor helicity method to give a concrete 
expression for the antenna function. Given values for the momenta ki, these quantities can be 
evaluated numerically. For example, 

K 2 



Ant(a",&- * — a+,l + ,6+) 
where I have used identities such as 



[a b] {ba) {a b) 
(al) (lb) (ab) 2 [ba] [ab] ' 



(3.6) 



(qa) = (ba) (qb) (a a) 
(qa) (ba) (qa) (ba) ' 
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(3.7) 



and have dropped terms such as the second which are non-singular in the limit. 

Squaring the antenna function, summing over final helicities (a, 1, b), and averaging over initial 



where the omitted terms are less singular, and in particular lead to finite integrals (as e — > 0) when 
integrated over 'singular' (soft or collinear) regions of phase space. (The inclusion or exclusion of 
such terms is arbitrary, as they are not universal. In practical applications one may choose to omit 
them.) The reader may verify that this expression reproduces known expressions in the soft and 
collinear limits. 

Away from four dimensions, we must take the unobserved momenta (a, 1,6) to have 4 — 2e 
rather than four components; depending on the variant of dimensional regularization, we may also 
take polarization vectors or observed momenta (a, b) to have 4 — 2e rather than four components. In 
the CDR scheme, one obtains the same expression (3.8) for the antenna function; in other schemes, 
there will in general be additional contributions. 

If instead we choose the Catani-Seymour reconstruction function, we cannot expect to repro- 
duce the correct behavior in the k b \\ k\ limit, but we will reproduce the Catani-Seymour dipole 
factorization in both the k a \\ k\ and k\ — ► limits. (Although the antenna function is not equal 
to the Catani-Seymour dipole function, the difference between them consists of terms which are 
not singular in the limit. Such terms, as noted above, are not universal, and may be omitted at 
will.) The corresponding helicity amplitudes could be thought of as the 'square root' of the dipole 
function. The antenna function may be thought of as a generalization of the dipole factorization 
formula. 

Although I have discussed only the antenna factorization function for purely gluonic ampli- 
tudes, it is clear that one can write down analogous formulae for mixed quark-gluon amplitudes. 
The approach presented here has a straightforward generalization to multiple-singular emission, 
which I shall endeavor to discuss elsewhere. I thank Z. Bern, S. Catani, L. Dixon, and N. Glover 
for helpful comments. 



helicities (a, b), one obtains 




(K 2 - Sab ) + sl h ] 2 



+ ••• , 



(3.8) 



S a lS lb S ab (K 2 ) 2 
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